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Particle production in deep inelastic scattering on nuclei is reduced due to absorption 
of the produced particles in the nucleus. The photon ejects a quark from a bound nucleon 
which propagates through the nucleus forming a prehadron before turning into a hadron. 
We calculate the higher twist effect in hadronization which dominates the z — > 1 region 
of fragmentation. 

1. Introduction 

Physics with nuclei extends the realm of strong interactions from the hadronic size scale 
of 1 fm to the size scale of nuclei of 10 fm. When a high energy probe illuminates the 
nucleus, a picture of the nucleus develops which is different from nuclear spectroscopy 
or low energy nuclear physics. Because of the high energy probe, QCD as the theory 
of quarks and gluons is needed. On the other side it is necessary to understand strong 
interactions at large distances. Three main issues are associated with hadronization in 
nuclei: 

• At which distances does color neutralize in the vacuum? 

• How can one see effects of prehadron formation? 

• Does the prehadron have a smaller cross section due to color transparency ? 

Renewed interest arises from recent studies of deep inelastic electron scattering on 
nuclei, where the produced hadrons are analyzed. The HERMES [QJI21II3I1] experiment 
has brought new insight into the question, since it has access to lower energy transfers 
v, where the hadronization occurs inside the nucleus, and to a larger range of fractional 
momenta z which determine the formation length of the hadron. The nucleus helps to 
track the space time evolution of a parton, since the nucleons play the role of very nearby 
detectors of the propagating object. Deep inelastic electron nucleus scattering has the 
additional advantage that the electron gives a well defined energy v to the struck quark 
propagating through cold nuclear matter. The understanding of this process is crucial 
for the interpretation of ultra-relativistic proton-nucleus and nucleus-nucleus collisions. 
Due to factorization in deep inelastic scattering the semi-inclusive cross section can be 
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Figure 1. Semi- inclusive hadron production in deep inelastic scattering on a target T in the 
pQCD factorization approach. Parton distribution functions (PDF) and fragmentation 
functions (FF) represent the non-perturbative input. 

described by the product of a parton distribution function (PDF) with a fragmentation 
function (FF) cf. fig. 1. In the parton model the probability that a quark with flavor / 
and momentum fraction x is present in the target is multiplied with the probability that 
it hadronizes into a definite hadron which carries a momentum fraction z of the quark, 
fig. 1 shows a schematic diagram of semi-inclusive deep inelastic lepton scattering (SIDIS) 
on a target T, and the definitions of the four momenta of the particles involved in the 
process. In SIDIS besides the scattered lepton I ' the leading hadron h formed from the 
struck quark is detected with energy Eh = zv in the target rest frame. The summation 
over flavors includes the product of the fragmentation functions and structure functions 
for each flavor. The experimental data on nuclear effects in hadron production are usually 
presented in terms of multiplicity ratios as functions of z,v or recently Q 2 : 



In the above definitions N A is the number of outgoing leptons in DIS processes on a nuclear 
target of atomic number A, while dN\/dz is the ^-distribution of produced hadrons in 
the same processes; the subscript D refers to the same quantities when the target is a 
deuteron. In absence of nuclear effects the ratio R\j would be equal to 1. In this paper we 
review the current status of the absorption model [ElE] in describing hadron production 
in deep inelastic scattering and add recent results on higher twist effects to semi-inclusive 
particle production at large z — > 1. 

2. Hadron Multiplicity Distribution 

The theoretical calculation takes into account several nuclear effects. Parton distri- 
bution functions and fragmentation functions both depend on the virtuality Q 2 of the 
DIS process. This adjustment to the scale Q 2 takes into account all radiated gluons 
before and after the photon quark interaction in the leading logarithm approximation. 




(1) 
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In nuclei gluon radiation may be affected by the partial deconfinement of color which 
follows from overlapping nucleons. Therefore in DGLAP evolution of nuclear structure 
and fragmentation functions the starting scale is smaller and they evolve over a larger 
interval in momentum compared with the corresponding functions in the nucleon at the 
same scale Q. Once color neutralization has taken place, gluon radiation stops as a source 
of energy loss. A color neutral prehadron has been produced which interacts with the 
color neutral other nucleons in the nucleus. In the work by A. Accardi, V. Muccifora, 
D. Griinewald and myself [ EJ El similarly to the gluon bremsstrahlung model [ [7j the 
time scale for prehadron formation is short for large z and small z. It increases with z 
because of Lorentz time dilatation and it is small for large z, since the hadron has to 
be formed instantaneously, otherwise the energy loss downstream is too large. In fig. 2 
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Figure 2. Schematic space time picture of hadronization in the LUND model. 

we show a schematic view of the hadronization process in the rest system of the nucleus. 
Our theoretical model follows closely the Lund model for the fragmentation process [|H]. 
The nucleons iVijJVjj, ... are taken at rest in the nucleus. They are the target particles 
for the developing fragmentation cascade of the struck quark. We differentiate between 
the prehadron formation points Pi and the hadron formation points F,. The space-time 
development of the fragmentation process begins when the quark q is ejected from a nu- 
cleon. The quark propagates and the colour string between the quark and the remnant 
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breaks into smaller pieces. 

Hadrons are ordered according to their rank i. Note that the first-rank hadron is always 
created at the end after the original quark has traversed a distance L 



where k = 1 GeV/fm is the string tension. This length L can be very large cf. fig. 2, but 
this does not mean that long strings exist in the nucleus, because the string breaks and 
prehadrons are formed. There are two relevant lengths for the fragmentation process, the 
average position (/*) at which the prehadron is formed and the average distance (lh) < L 
at which the hadron is formed, where the averaging is over all hadron ranks. In the Lund 
model the formation length of the prehadron and hadron are related [|S]: 



They both increase linearly with the virtual photon energy v. However, as functions of z 
they behave rather differently, especially at z — > 1, where (/*) — > and (lh) — > L. 
Once the prehadron is produced, it makes inelastic collisions which degrade its z fraction. 
Since the fragmentation function is very steep, we consider the prehadron lost after an 
inelastic collision. We introduce a survival probability Sa(z, v) which is calculated in linear 
transport equations which include formation and absorption of the quark and prehadron. 
These equations improve previous work by Bialas and Chmaj [E] by treating hadronization 
and absorption on equal footing. The survival probability represents the probability that 
the prehadron and the hadron survives the traversal of the nucleus. 
We calculate the multiplicity distribution of a hadron h with momentum fraction z in a 
nucleus A with the following expression 



the differential cross section for a 7*g scattering computed in pQCD at leading order. 
The nuclear structure function £,a(Q)Q 2 ) and the nuclear fragmentation function 
D 1 j{z 1 ^a{Q)Q 2 ) are used. In ref. [ |H] we fitted a prehadron cross section equal to the 
2/3 of the hadron cross section which is realistic since the prehadron coming from gluon 
radiation of the struck quark is smaller than the final hadron. 

3. Prehadron absorption 

The dependence of hadron attenuation on the mass number A of the target nucleus is 
commonly believed to clearly distinguish the absorption and energy loss mechanisms for 
hadron attenuation. Indeed, the common expectation is that attenuation in absorption 
models is proportional to the hadron in- medium path length d, leading to Rm oc A 1 / 3 . 
On the other hand, the average energy loss for a parton traversing a QCD medium is 
proportional to d 2 , which leads to Rm oc A 2 ! 3 . 



(k) = (I*) + zL . 



(3) 
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Figure 3. Distribution of prehadrons and hadrons with z = 0.1, z = 0.5, z = 0.9 in Kr 
induced by an interaction of a photon with a corresponding average energy (u)(z) given 
by the HERMES [UJI21E1II! experiment. 



In fig. 3, we show a picture where the different formation areas of a premeson (left part 
of the figure) and a meson (right part of the figure) are shown relative to a target Kr 
nucleus which has a radius of about 5 fermi. From top to bottom the figures represent 
different z- values, z = 0.1, z = 0.5 and z = 0.9 are shown. The energy of the initial quark 
for the different values of z is given by the respective average of the virtual photon energy 
(v){z) given by the HERMES [UJ 121 El H] experiment. One sees that the hadron dom- 
inantly appears outside of the nucleus, whereas the prehadron forms inside the nucleus, 
especially for large and small z-values. Therefore we take the model described in ref. [ 
ITU] in the one step (= prehadron) approximation and neglect final hadron production 
in the calculation of the survival probability. To further simplify, we consider the case 
of a hard-sphere nucleus of mass number A and radius R = r^A 1 ^ 3 , with r$ = 1.12 fm. 
Neglecting absorption in deuterium and the small rescaling correction we find that the 
hadron multiplicity ratio Rm equals the hadron survival probability Sa- These simplifi- 
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cations make an analytical computation of the attenuation possible (cf. ref. [El) and an 
expansion of the result in terms of the mass number yields 

1-Rm = ciA 2/3 + c 2 A 4/s + 0[A 2 ] . (5) 

The numerical values of Ciand c 2 depend on z and u, in general c 2 « C\. A fit to 
the theoretical model without any simplifications and to the experimentally measured 
multiplicity ratios finds 1 — R m = cA a with a strong correlation between a and c favoring 
a = 2/3 in agreement with the above considerations on the absorption model. 

4. Higher Twist Effect near z — > 1 

In the late eighties E.L. Berger [|TTj has published several papers on higher twist effects 
in semi-inclusive particle production. The fragmentation function Dh/ q (z, Q) contains a 
leading twist term D^J q (z,Q) which evolves logarithmically with Q 2 due to perturbative 
gluon radiation and a higher twist term D^(z, Q) responsible for direct hadron produc- 
tion 

D h/q (z,Q 2 ) = Dff q (z,Q) + Dff q (z,Q). (6) 

The higher twist contribution to the fragmentation function for production of pions is 
modeled in ref. [E] by: 

D 7r/q (z,Q 2 ) = D^(z,Q)+ 1 -F 7T (Q 2 ). (7) 

where F n (Q 2 ) is the pion monopole form factor F n (Q 2 ) oc m 2 /(m 2 + Q 2 ) determined 
by the p meson mass. The leading twist fragmentation function at low Qq = 2 GeV 2 
behaves as D^J q (z, Qq) oc (1 — z). This means that for z > 0.7 the higher twist direct pion 
production takes over. The radiated gluon and remaining quark must carry each about 
one half of the momentum of the struck quark. The resulting propagators suppress the 
higher twist contribution. 

For nucleon production, we expect an even larger higher twist suppression factor: 
D N/q (z,Q 2 ) = D$ q ( z ,Q)+D%J q {z,Q) 

= D L J /q {z,Q) + l -F N (Q 2 ). (8) 

The resulting nucleon form factor has a dipole behaviour. A comparison of our previous 
calculations [ 6 with the experimental data indicates that for large z- values an additional 
mechanism is at work. Therefore we have estimated that the attenuation of the directly 
produced higher twist hadron must be less , due to color transparency. We fitted a reduced 
cross section a CT = 0.37a 7rN = 10 mb both for the direct pion and nucleon production. In 
fig. E|we show the results for the multiplicity ratios Rm- In both cases Rm falls off more 
slowly for z — > 1. For the nucleon it more or less flattens out. 

This purely phenomenological study puts numerical estimates by Kopeliovich and Nemcik 
into a new light. Indeed for z > 0.7 direct hadron formation is important for small Q 2 . 
Note the average Q 2 = 2.5 GeV 2 in the Hermes data. The formalism to be developed 
may allow to calculate the color transparency effect both for the leading twist prehadron 
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Figure 4. Multiplicity ratios Rm{Kt) for pions (left panel) and protons (right panel) as 
a function of z. Dashed lines show the leading twist fragmentation process of ref. [ |H| 
computed in the one step approximation. The full lines show the modifications due to 
the higher twist process. The points show the experimental HERMES data 



and the higher twist produced direct hadron. Deliberately we left out to parametrize the 
Q 2 dependence of the color transparency reduced cross section a CT , since the data on 
Rm do not need such a dependence. On the other side the semi-inclusive cross sections 
themselves should show a strongly Q 2 dependent higher twist effect in the region z > 0.7 
and Q 2 < 4GeV 2 . Higher twist fragmentation is additionally enhanced in the nucleus 
because of its smaller absorption cross section and thereby becomes better accessible 
experimentally. 



5. Conclusions 

In summary, the absorption model [03 E] can describe the HERMES data on the nuclear 
modification of hadron production in DIS processes. A reasonable choice of parameters 
with the vacuum string tension and a prehadronic cross section equal to 2/3 of the hadronic 
cross section explains the data. The A-dependence in the absorption model is simple. The 
attenuation 1 — Rm is proportional to A 2 / 3 since mostly prehadrons from the away surface 
of the nucleus contribute to the cross section. 

Finally high energy electron-nucleus collisions allow to observe the higher twist effect in 
the fragmentation function at z — > 1. The directly produced higher twist hadrons have 
a smaller prehadron cross section a = 10 mb due to color transparency. The nucleus 
enhances this higher twist effect. 
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